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ABSTRACT 


In a quest to improve space-based observational ca- 
pability, an increasing number of investigators are 
proposing missions with precision formation flying 
architectures. Typical missions include the Micro- 
Arcsecond X-ray Imaging Mission (MAXIM), Stel- 
lar Imager (SI), and the New Worlds Observer 
(NWO). Missions designed to explore targets in deep- 
space generally require holding a formation config- 
uration fixed in inertial space during science obser- 
vation. Analysis in this paper is specifically aimed 
at the NWO architecture, characterizing the natu- 
ral drift of the line-of-sight and the separation range 
for two spacecraft operating in the vicinity of the 
Earth/Moon-Sun L 2 libration point. Analysis em- 
ploys a linear form of the relative dynamics associated 
with an n-body gravity field. The study is designed to 
identify favorable observation directions, character- 
ized by minimal line-of-sight drift, along the mission 
timeline. 


NOMENCLATURE 


A(t) = 

Dynamics matrix for linear form of 
equations of motion 

r. 

Position vectors, subscripts depicted 
in figures 

U thruat,F ~ 

External control force applied to 
Follower spacecraft 

'^'thrxiatjL 

External control force applied to 
Leader spacecraft 

X = 

Position of Follower referenced to 
Leader position 

I 

V = 

Superscript designating inertial (I) 
frame 

Pi “ 

Gravitational parameter for i th body 
of n-body system 

11*11 

The 2-norm of the vector x 


1. INTRODUCTION 


Precision formation flying characterizes missions requir- 
ing relative state maintenance to very tight tolerances. 
Consequently, precision formation flying missions typi- 
cally include continuous control schemes. Many candi- 
date precision formation flying architectures are derived 
for missions designed to explore targets in deep space. 
These missions may be characterized by inertially fixed 
formation configurations during science observation pe- 
riods. Generally, control tolerances are relaxed during re- 
configuration maneuvers. Typical missions include: the 
Micro-Arcsecond X-ray Imaging Mission (MAXIM) [1], 
Stellar Imager (SI) [2], and the New Worlds Observer 
(NWO) [3]. 

NWO is a survey mission designed to discover and ex- 
plore potential earth-like planets that exist outside the so- 
lar system. The mission concept includes two spacecraft, 
a telescope and an occulter. The telescope is designed 
to observe a planetary system about a star. The occul- 
ter is designed to block light from the star allowing ob- 
servation of much dimmer planetary bodies orbiting the 
star. The challenge is to maintain the occulter on the 
telescope/target star line-of-sight, and to a lesser degree 
maintain the required telescope/occulter separation. With 
this configuration performance of NWO is much more 
sensitive to drift of the occulter off the line-of-sight di- 
rection than drift along the inter-spacecraft range. 

To elaborate, assume the telescope is properly aimed 
along the line of sight to the target star. Optical perfor- 
mance requires the telescope/occulter spacecraft line of 
sight be maintained within an angle of 5 milli-arcseconds 
at 70,000 kilometers, equivalent to maintaining the oc- 
culter position within a meter (approximate) of the tele- 
scope/target line of sight. In contrast allowable drift 
along the line of sight is on the order of 100 kilome- 
ters from the desired range. Clearly performance of the 
telescope/occulter is much more sensitive to drift normal 
to the telescope/target line of sight than along the line 
of sight, relative spacecraft range. This motivates a de- 
sire to align the natural drift between the two spacecraft 
along the separation vector, minimizing the line-of-sight 
drift. Previewing this paper, the approach exploits the lin- 
ear form of the relative dynamics, identifying the eigen- 




vectors for the dynamics matrix in the linear state equa- 
tion. Minimal drift in the direction normal to the tele- 
scope/target line of sight is achieved by aligning the ob- 
servation vector with an eigenvector of the matrix repre- 
senting the linear form of the relative spacecraft dynam- 
ics. The eigenvectors evolve with the formation trajec- 
tory. Ideally with a target rich environment, a targeting 
strategy might constrain observations to lie along (near) 
a particular eigenvector. In practice this strategy severely 
limits mission planning and target selection. A better 
strategy weights relative spacecraft drift as one of the fac- 
tors in designing the mission observation schedule. 

True coalignment of an eigenvector and the line-of-sight 
is not achievable for the duration of a given observation. 
The linear form of the dynamics are time varying. The 
eigenvectors of the dynamics matrix slowly rotate in in- 
ertial space. The observation vector is fixed in inertial 
space. 

While the analysis centers on the NWO mission, the de- 
velopment and results generically apply to any similar 
formation architecture. 


2. RELATIVE DYNAMICS 


As noted, the analysis goal is to identify the relative range 
vector orientations that minimize natural drift of the 
range vector direction. The natural relative dynamics for 
a spacecraft formation operating near the Earth/Moon- 
Sun L 2 libration point are principally forced by differ- 
ential gravity and solar pressure. The gravitational differ- 
ential acceleration is a space environment characteristic, 
independent of spacecraft design. In contrast, differen- 
tial acceleration due to solar pressure is a design consid- 
eration tightly correlated to the physical properties of the 
spacecraft, and subject to variation with the orientation of 
the spacecraft with respect to the incident solar radiation. 
This analysis is limited to characterizing the differential 
gravity field. 

Reference [4] presents the development of the linear form 
for the relative dynamics in inertial coordinates for the 
restricted three body problem, depicted in Fig. 1. Refer- 
ences [5] and [6] extend the development for the n-body 
problem, depicted in Fig. 2. The development of the lin- 
ear form of the relative dynamics for the n-body problem 
is repeated in summary form for reference. 


2.1. Dynamics of Relative Motion for the n-Body 
Problem 



Figure 1. Two Spacecraft Orbiting in the Earth/Moon - Sun 
Rotating Frame 



Figure 2. Two Spacecraft Under the Gravitational Influence of 
n-Bodies 


Treating each body as a point mass and allowing space- 
craft thrust, the equation of relative motion of the Fol- 
lower with respect to the Leader is: 



The n-body problem examines the behavior of an in- 
finitesimal mass in the combined gravitational field of ’n’ 
finite masses orbiting their common center of mass. A 
typical two spacecraft formation operates in the gravita- 
tional field of n-bodies, as shown in Fig. 2. 


Eq. 1, presented for a gravity field in an inertial frame, 
provides an exact expression of the nonlinear dynamics 
of relative motion between the Follower and Leader 
spacecraft. The next step expresses the relative dynamics 



of the Follower with respect to the Leader in linear form. 


2.2. Linear Form of Relative Dynamics 

Transforming Eq. 1 to a linear form requires examination 
of each term. Consider j ~ } 
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Apply a binomial expansion to first order, resulting in 
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It follows from Eq. 3 that 
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Replace ( x T x ) r iL and (r iL J x) r iL in Eq. 5 with the 
equivalent expressions, (r iL x T ) x and (r iL r iL J ) x, 
yielding 
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In order to remove the state dependence of the 
dynamics matrix in Eq. 6, assume a tight for- 
mation with ||x|| << ||r iIf ||, for all i. Thus, 

<< 1. Likewise, the 

11 ’ iL ii || 

term | ^j 2 is also neglected, since its matrix norm is 

much smaller than the remaining terms. Eq. 6 simplifies 
to 
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Note that denotes a unit vector along r iL 


Substituting Eqs. 3 and 4 into Eq. 1 gives 


In summary the linear form of the dynamics is expressed 
as 



= E(t) x + u t 


( 8 ) 


where 

's(t) =-£p^[l3 + 3(e 4i e it T )] 


The dynamics in matrix form are 


£ = 1 A(t) £ +B (u th _ t>F -u thra . t ' L ) (9) 


Where e SL (to) and e EL (to) represent unit vectors from 
the Sun to the Leader, and the Earth/Moon mass center 
to the Leader, respectively. By inspection one eigenvec- 
tor is a unit vector in the direction of the cross product of 
e SL (t o) and e EL (t 0 ) with an eigenvalue of — (ci 4- c 2 ). 
The two remaining eigenvectors lie in the plane spanned 
by e SL (t 0 ) and e £L (t 0 ). The special case with e SL (t 0 ) 
and e EL (to) collinear is not considered since the config- 
uration is not achievable along expected trajectories for a 
mission near the Earth/Moon-Sun L 2 point. The typical 
orientation of the eigenvectors is depicted in Fig. 3. For 
an orbit about the Earth/Moon-Sun L 2 point the eigen- 
value associated with e 2 is positive, representing an un- 
stable mode. The eigenvalues of e 1 and e 3 are negative, 
representing stable modes. Stability properties are for the 
nominal S (t 0 ). 


Where: 






I 3 

0 


B = 


0 

I 3 


2.3. Characterizing Relative Drift 


Based on the previous discussion, the term S (t) char- 
acterizes the gravity gradient between the two spacecraft. 

S (t) provides a linear form of the relative acceleration, 

referenced to the leader spacecraft. S (t) is time varying 
with a very slow evolution correlated with the nominal 
6-month period for an orbit about the Earth/Moon-Sun 
L 2 point. Using the NWO mission concept, the science 
schedule generally limits the total time for each individ- 
ual target to two days with less than a day of actual sci- 
ence observation time. Consequently, the relative space- 
craft dynamics for a single observation can be character- 
ized by a nominal, constant S (t). 

Following the assumption *3 (t) is constant, the magni- 
tude and direction of the relative spacecraft acceleration 
is defined by the eigenvalues and associated eigenvectors 
of 3 (t 0 ). t 0 represents the nominal point in time for the 

science observation. Referring to Eq. 8, 3 (t) is sym- 
metric. Therefore, the eigenvectors are mutually perpen- 
dicular. The orientation of the eigenvectors is approxi- 
mated by evaluating 3 (t 0 ) with only the contributions 
of Sun and Earth/Moon system gravity fields, expressed 
as: 



Figure 3. Orientation of Eigenvectors (e 1 , e 2 , e 3 J of 
H(t 0 ) relative to e SL (t 0 ) and e EL ( t 0 ) 


Alignment of the formation along any of the eigenvector 
directions restricts the relative spacecraft (gravitational) 
acceleration along the separation vector. Following the 
eigenvector designations from Fig. 3, the e x direction 
generally aligns nearly with the spacecraft to Sun line, a 
direction typically prohibited for science observations by 
deep-space exploration mission constraints. Therefore, 
e 2 and e 3 are the allowable, principal observing direc- 
tions that minimize angular drift of the line-of-sight. 


3. SIMULATION 
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Simulation results demonstrate the variation in the radial 
and crosstrack line of sight drift along a representative 
trajectory for the NWO mission. The observation vec- 
tor is fixed in inertial space with a spacecraft separation 
range of 72,000 kilometers. The nominal orbit for the 
telescope spacecraft about the Earth/Moon-Sun L 2 point 
is shown in Fig. 4. 


ci 
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Figure 4. Leader (Telescope) Spacecraft Trajectory about 
Earth/Moon-Sun L 2 point 


3.1. Scenario 1 - Observation Vector Normal to 
Ecliptic Plane 

In the first scenario the observation vector is defined to 
align with the normal direction for the ecliptic plane. 
Fig. 5 shows the relative acceleration between the two 
spacecraft with respect to the mission timeline. The 
principal component of the relative acceleration remains 
aligned with the observation (range) vector along the en- 
tire orbit. 



Figure 5. Relative Acceleration of Follower with Re- 
spect to Leader for Observation Vector Normal to Eclip- 
tic Plane 


3.2. Scenario 2 - Observation Vector within Ecliptic 
Plane 

In the second scenario the observation vector is defined 
to lie within the ecliptic plane. The observation vector is 
initially aligned with the Earth/Sun line. Fig. 6 shows the 
relative acceleration between the two spacecraft with re- 
spect to the mission timeline. The components of the rel- 
ative acceleration exhibit significant variation along the 
mission timeline. Fig. 7 depicts the same results as a 
function of the angle between the observation vector and 
the eigenvector for the unstable mode of S (t). The un- 
stable mode is the e 2 direction shown in Fig. 3. The 


results indicate the crosstrack drift is minimum for ob- 
servation vectors along the unstable mode eigenvector of 
E(t 0 ), or perpendicular to that vector. The results apply 
to a vector confined to the ecliptic plane. NWO mission 
science constraints prohibit observations along the unsta- 
ble mode eigenvector, e 2 , since it nominally aligns with 
the spacecraft to Sun vector. 



Figure 6. Relative Acceleration of Follower with Re- 
spect to Leader for Observation Vector within the Ecliptic 
Plane as a Function of Mission Timeline 



Figure 7. Relative Acceleration of Follower with Re- 
spect to Leader for Observation Vector within the Ecliptic 
Plane as a Function of Angular Separation Between Ob- 
servation Vector and eigenvector e 1 of S (t 0 ). 


4. CONCLUSIONS 


This work demonstrates the utility of employing the lin- 
earized relative dynamics for a spacecraft formation as 
a planning tool for scheduling observations. The lin- 
earized dynamic model identifies preferred observation 
directions which minimize natural drift of the formation 
line-of-sight. Preferred directions align with the eigen- 
vectors of the matrix representation of the linear form of 
the dynamics. The strategy assumes the mission science 
goals are achievable with the imposed constraint on the 
observation line-of-sight. 
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